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r-| . Abstract 

k> ' ^ e sri ow that in the minimal three generation seesaw models for neutrinos, the presence of leptonic 
^ . (L e + L^ — L T ) x S2 symmetry leads to one of the right handed neutrinos remaining massless. This 
5^ ! state can then be identified with the sterile neutrino required for a simultaneous understanding 
of solar, atmospheric and LSND observations. We present a gauge model where the presence of 
higher dimensional operators originating from Planck scale physics lead to a realistic 2+2 mixed 
scenario that fits all oscillation data. The model predicts a range for the mixing angle U e % and an 
effective mass for neutrinos emitted in tritium decay, which can be used to test this model. 
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1 Introduction 



The evidence for neutrino masses and mixings now appears to be quite solid from the positive 
results for neutrino oscillations from solar neutrino data in seven different experiments, Chlorine, 
Kamiokande, Super-Kamiokande, SAGE, GALLEX, GNO, SNO|I| and atmospheric neutrino data 
from Super-Kamiokande, 1MB, Soudan and MACRO experiments ||. The third piece of evidence 
is from the Los Alamos LSND experiment that shows an oscillation from the muon neutrino to 
the electron neutrino (from z^j to Pe)||. The KARMENQ experiment which looked for 17^ to TF e 
oscillation did not find any evidence for it and eliminated a large fraction of the parameter space 
allowed by LSND. It is hoped that the Mini-BOONE experiment at FERMILABg] will settle the 
issue in near future. 

Theoretical analyses have made it clear that a simultaneous explanation of all the neutrino 
oscillation observations requires the existence of an ultralight sterile neutrino in addition to the 
three known active neutrinos (z/ e , u^, z+)0- There are two mass patterns that are currently under 
active discussion; the 2+2 schemeQ and 3+1 scheme @. The 2+2 scheme has two mass eigenstate 
neutrinos with mass around an eV and a lighter pair with mass near 10~ 3 eV, with the latter pair 
explaining the solar neutrino data, the former explaining the atmospheric neutrino data and the 
gap between the two pairs explaining the LSND results. Recent SNO data disfavors the original 
version of the 2+2 model where all the missing solar z/ e 's are converted via a small angle MSW 
mechanism only to the sterile neutrinos. The situation where only a fraction of the missing v e s 
convert to u s s and the rest to active ones has been studied in detail|| and shown to be quite 
viable at the moment. The 3+1 picture on the other hand is severely constrained by the known 
accelerator data but is also viable for certain values for Am 2 . 

If indeed the sterile neutrino turns out to be required, an important theoretical challenge is to 
seek an understanding of such tiny mass of a particle which is neutral under the standard model 
forces. This is in clear distinction from the case for the active neutrinos for whom the small mass 
has a very convincing explanation in terms of the seesaw mechanism^ that, at the basic level 
requires only that three singlet right handed neutrinos be introduced with masses close to the GUT 
scale. The conventional seesaw mechanism makes the three right handed neutrinos superheavy, 
leaving the three active neutrinos with a very light mass, as seems to be case. In addition to 
this appealing feature, this minimal seesaw model also restores quark lepton symmetry into the 
fundamental fermions of nature, making this the standard paradigm of neutrino mass physics. 
Generic attempts to understand the ultralight sterile neutrino within the seesaw framework involve 
the introduction of a new singlet fermion into the theory beyond the three right handed neutrinos 
just discussed and use extra symmetries to protect it from being superheavy. While there exist 
perfectly acceptable, interesting and technically natural scenarios for the sterile neutrino with 
extra singlet fermions|10[, these extra assumptions have caused skepticism regarding the existence 
of the sterile neutrino and in turn about the LSND results. 

An interesting theoretical challenge is therefore to explore whether within the minimal quark 
lepton symmetric seesaw framework, one can have an ultralight sterile neutrino without the need 
for additional singlet neutrinos, outside the quark-lepton symmetric framework and fit all obser- 
vations including LSND. Such a project is apriori not implausible in view of the fact that the right 
handed neutrinos already present in the quark-lepton symmetric framework have all the proper- 
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Table 1: Relevent right handed fermion and scalar fields and their transformation properties. Here 
we have defined Y = I 3 r + 



ties of a sterile neutrino, except the mass. In a recent paper |TT| it was shown that if there exist 
certain leptonic symmetries in the seesaw model, they can make one of the right handed neutrinos 
massless in the tree approximation of a renormalizable theory making it possible to identify it 
with the needed ultralight sterile neutrino. Using this idea for the case of an (L e — — L T ) x 5 Mr 



symmetry, (where S^ T is the interchange symmetry between and Vtr), it was shown in |TT 
that the 3+1 scenario for the LSND observations can be reproduced within the minimal seesaw 
picture. 

In this letter, we show that when the leptonic symmetry is instead chosen to be (L e + L M — 
L T ) x S2, the suggestion of ref.|ll[] can be extended to obtain the 2+2 explanation of LSND with 
one of the right handed neutrinos of the seesaw picture playing the role of the sterile neutrino. The 
mixing of the sterile neutrino to the active neutrinos arises from Planck scale suppressed higher 
dimensional operators. A consistent picture with 10 13 GeV for the seesaw scale and 10 6 GeV for 
the breaking of the leptonic symmetries seems to emerge and leads to mass matrix ansatz for the 
2+2 case recently proposed in JL2J. Experimental tests of this idea are proposed. 



2 Seesaw with a leptonic symmetry and right handed neu- 
trino as an ultralight sterile neutrino 

In this section, we discuss how imposing (L e + — L T ) x S eft in a seesaw framework leads to a 
massless right handed neutrino which can play the role of the sterile neutrino. We will work with 
the gauge group SU(2)l x U(1)i 3R x U(1)b-l with particle assignment as in Table I. We call this 
the minimal seesaw picture since all the fermions in our model are needed in the conventional 
seesaw explanation of the small masses of the neutrinos. The Higgs field A transforming as 
(1,-1, +2) under the gauge group couples to the right handed neutrinos (f e jj, v^r, v t r) with 
couplings typically of the form z/ri/rA. After spontaneous symmetry breaking, the field A acquires 
a vev i.e. < A >= vr and breaks the gauge symmetry to the standard model. This gives mass to 
the right handed neutrinos and is the seesaw scale. The Higgs field 0(2, 1/2, 0), whose weak scale 
vev breaks the standard model gauge group, gives Dirac mass to the neutrinos. If we denote the 
mass matrix for the right handed neutrinos as Mr, then the complete 6x6 mass matrix involving 



3 



the Dirac mass and Majorana mass for the neutrinos can be written as 
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When none of the eigenvalues of Mr vanishes, one can obtain the mass matrix for the light 
neutrinos as 



mIMr x m, 
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(2) 



This is the so-called type I seesaw formula. On the other hand when Mr and Mr, matrices have 
zero eigenvalues, one must "take them out" of the matrix before using the seesaw formula. As was 



noted in [pH| , this turns out to be the case when there are leptonic symmetries. To obtain the 2+2 
scenario, we will consider the leptonic symmetry to be (L e + L M — L T ) x S^r. At the beginning of 
our analysis we go to a basis where charged lepton mass matrix is diagonal where leptonic mixing 
matrix is simply neutrino mixing matrix. In this basis we identify the leptonic flavor indices e, /x, r. 
After that the renormalizable Lagrangian invariant under (L e + L M — L T ) x S^r symmetry has the 
form 
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Spontaneous symmetry breaking then leads to Dirac and Majorana type mass matrices of the 
form, 
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The 6x6 mass matrixflll would then take the form 
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In the rotated basis it is 
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where 
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k' + k 2 
kv e L + k'vl 



k' 2 + k 2 
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(6) 



(7) 



(9) 



and £ = 2 ^ kk ' , We next employ the see-saw mechanism to generate the light Majorana neutrino 
mass matrix, 
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In the limit k = k! we have £ = 2k and m = v ^^ n33 |TT[| , however because S 2 l symmetry is absent 
we have the freedom to choose k 7^ k' which will help to achieve small mixing in the u e — sector 
compatible to LSND. 

The challenge after this is to induce large left right mixing. There are higher dimensional oper- 
ators in this scenario which are sufficient to induce left-right mixing after the see-saw mechanism 
as will be described below. We can write down following higher dimensional operators 



L e v„ R (a )(4>) 
L^v- R (<r )(ij>) 



(11) 
(12) 

(13) 

(14) 
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Now the 4x4 mass matrix in low energy can be written in terms of rotated low energy fields. 
We identify v% = u' e , v L = v'^ v_ R = v s 
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CASE 1 



We notice that in the limit m — > 0, 5 — > the eigenvalues are 
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Thus LSND mass splitting squared can be of order M 2 . Introduction of non-zero m and 5 will 

l atm- To leading order the eigenvectors are, 
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Depending on the hierarchies between mi,m2,m 3 ,5 and m various scenarios can emerge from 
general expressions given in Eqn. fllTl). Let us consider a simple possibility. Taking mi » 
mi, 777.3, m i an d mi ~ 7ri 3 we have M = m 2 . Eigenvectors reduce to 
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This means that the massive pair is constituted mainly of v' and u s , whereas the light pair is 
constituted mainly of u' e and v T and each pair has large mixing among them. Now take k >> k' 
in Eqn. (|7|) and Eqn. (|3|). Then z/ « z/£ and z/ ~ z/£. We must introduce non-zero m, 5 mi and 
7?72 in such a way that the two pairs mix among each other to explain the LSND result. 

This can be done numerically for example if we take the parameter set, m = 0.04, mi = 



0.05,m 2 = l,m 3 = 0.06,(5 = 0.02 eVs. We get, Am 2 & = 2 x 10" 5 , Am 2 atm = 4 x 10" 2 , Am 2 
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(19) 



We note that now we have introduced small v e — v„ mixing even though the primary mixing pattern 



is governed by Eqn. (18). We can now calculate P( 



for LSND experiment to be 



>LSND 



On the other hand, we have, 
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Combining these two equations we have 



sin 2 26 LSND = 0.01 

Thus we see that we get Am 2 as well as sin 2 29 compatible with LSND experiment] 



we have already got large (near-maximal) v e 



(22) 
. Whereas 

«-> v T and large (near-maximal) «-> v s mixing to 
account for solar and atmospheric neutrino deficits. However in two-generation analysis of the 
atmospheric neutrino data the pure «-> v s oscillation scenario is seen to be highly disfavored, 
as it fails to reproduce the correct zenith angle distribution of the data, the effect being more for 
the higher energy neutrinos. The sterile option is also in contradiction with the neutral current 



sample of the SK atmspheric neutrino data [FJ| . But from the four- neutrino global analysis of the 
solar and the atmospheric neutrino data taken together || this scenario is not ruled out, although 
disfavored ||. The other drawback of this mass texture is that it predicts near-maximal mixing 
for the solar v e <-> v T oscillations. However, after the inclusion of the SNO data, maximal mixing 



in the LMA region is seen to be disfavored in most analyses of the solar neutrino problem |14 



CASE 2 



Let us apply the inverse of the rotation given in Eqn (|7]) and Eqn (^) on the basis vectors 



[v e , u' ). Then matrix in Eqn. ([TJ) can be expressed as, 
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For small b this mass matrix has been studied in Ref. 

one discussed above. It was shown in detail that it is suitable for the mixed 2 + 2 scenario. 



T^| using a different symmetry than the 

A 



similar matrix was also obtained under radiative scheme. However the 44 diagonal element vanish 



genericaly in simple radiative schemes ||1 5|[ 
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3 Predictions of the Model 



With the range of values for the parameter set needed to explain the global solar, atmospheric and 
the accelerator /reactor data including LSND, our four-generation model makes definite predictions 
for the effective mass of the neutrinos in beta decay and neutrinoless double beta decay processes. 
For the neutrinoless double beta decay the "effective Majorana mass parameter" is given by 



< m > 



i=lA 



mJJ„, 



(24) 



the 11 element of the mass matrix. Since this element is identically zero for the model under 
consideration, our model predicts no neutrinoless double beta decay. In a recent paper [|1(J has 
claimed of a positive evidence for neutrinoless double beta decay, based on the reanalysis of 
the data from the Heidelberg-Moscow experiment. A best-fit Majorana mass of 0.39 eV for the 
neutrinos at 95 % confidence level was reported. However there are proposals for new generation 
neutrinoless double beta decay experiments which are expected to have better sensitivities so they 
can resolve this issue at higher confidence levels. These include NEM03 
GENIUS 



17 , CUORE [18 and 



B3 



which hopes to probe neutrino masses down to ~ 10~ 2 eV. This is a good way to 
falsify the scenario we are presenting here. 

The beta decay experiments look for evidence of neutrino mass by observing the end point 
spectrum of 3 H beta decay. The best bound available to date is pO, |21j 



m Ue < 2.2 eV at 95% C.L. 



(25) 



In terms of the mass and mixing angles the quantity probed by the beta decay experiments is 
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For the CASE 2 of our model we have to linear order in 5, e\ and e 2 
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With the range of values for the parameters Wq, a, €\ and e 2 required to satisfy the global neutrino 

m„„ ~ 0.03 for Am 2 LSND ~ 1 eV 2 . The forthcoming 

0.3 eV |22|]. Thus we expect no 



oscillation data we expect the value of 



KATRIN experiment is expected to have a sensitivity upto m Ue 
signal for m Ue in KATRIN from our model. This is another way to test our model. 

The other intriguing quantity that needs precise measurements is \U e s\ 2 + |£/ e 4| 2 which will give 
the fraction of v e mixed in the upper doublet. For the CASE 2 upto linear terms in 5, e± and e 2 
we have 
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while 

m 3 u, 3 + u* e4 u, 4 \ 2 « ^(tt^) 2 (29) 

« sin 2 29 LS nd 

Thus for the CASE 2 of our model |£7 e3 | 2 + |f/ e4 | 2 « (sin 2 20 LSN d)/± ~ 0.0007 much more stringent 
than the current bound from the short baseline terrestrial experiments (BUGEY[f23f) which is 

\U e3 \ 2 + \UeA? -0.01 (30) 

Values of |t/ e 3| 2 + |^e4| 2 as low as 0.0007 would be difficult to detect with usual laboratory beams. 
However with proposals of high intensity neutrino beams from neutrino factories, it would be 
interesting to see if one can probe to such accuracy. 



4 Comments and conclusion 

A few comments are now in order: the first theoretical comment is on the choice of the gauge 
group: while we have illustrated our basic idea in the context of an SU(2)l X U(l)i 3R x U{1)b-l 
gauge group, it could easily be implemented in the context of the standard model gauge group 
as well. However, in this case the meaning of the seesaw scale will remain a mystery. One could 



also use the left-right symmetric gauge group SU(2) L x SU(2) R x [/(I)b-lPH- In this case since 
the right handed neutrino and the right handed charged leptons are in the same multiplet, S efl 
symmetry leads to electron being mass less prior to symmetry breaking. One can however extend 
the model to make the electron massive. Quark lepton unification of these models is nontrivial 
due to absence of any analog of L e + L M — L T symmetry in the quark sector. 

The breaking of the leptonic symmetries can be made soft by adding to the Lagrangian terms 
such as (cr^~) 2 so that there are no massless scalar bosons. 

In conclusion, we have presented a minimal seesaw model using three active and three right 
handed neutrinos (as would be suggested in a quark- lepton symmetric theory), which can lead 
to an ultralight sterile neutrino. The sterile neutrino in this case is none other than the right 
handed neutrino. Spontaneous breaking of the leptonic symmetries lead to mixing of the sterile 
neutrino with the the active neutrinos resulting in a 2+2 hybrid scenario for LSND, which provides 
a viable description of neutrino oscillation data. The model predicts an effective U e3 ~ 0.025 and 
an effective end point mass in tritium beta decay ~ 0.03. 
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B.B would also like to thank Amitava Raychaudhuri and Probir Roy for discussions on Ref. [ 1 1 



References 

[1] Y. Suzuki et ai, Super-Kamiokande collaboration, Nucl. Phys. Proc. Suppl. 77, 35 (1999); 
B. Cleveland et ai, Ap. J. 496, 505 (1998);J. N. Abduratshitov et al., SAGE collaboration, 



9 



Phys. Rev. C 60, 055801 (1999); W. Hampel et al, GALLEX collaboration, Phys. Lett. 
B447, 127 (1999); M. Altman et al, GNO collaboration, Phys. Lett. B490, 16 (2000); Q. 
Ahmed et al, SNO collaboration, nucl-ex /0 1060 1 5| . 



[2] Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 (1998); R. Becker-Szendy et al., Phys. Rev. D 
46, 3720 (1992); S. H. Hirata et al., Phys. Lett. B 280, 146 (1992); W. W. M. Allison et al, 
Phys. Lett. B 449, 137 (1999); M. Ambrosio et al., Phys. Lett. B 517, 59 (2001). 

[3] C. Athanassopoulos et al, Phys. Rev. C54 (1996) 2685; C. Athanassopoulos et al.,, Phys. 
Rev. C58, 2489 (1998); A. Aguilar et al, |hep-ex/ 0104019. 



[4] K. Eitel et al., KARMEN collaboration, |hep-ex/ 0008002 



[5] See talk by J. Conrad, WIN2002 workshop in Christ Church, Newzealand (January, 2002). 

[6] D. O. Caldwell and R. N. Mohapatra, Phys. Rev. D 48, 3259 (1993); J. Peltoniemi and J. 
W. F. Valle, Nucl. Phys. B406, 409 (1993); J. Peltoniemi, D. Tomassini and J. W. F. Valle, 
Phys. Lett. B 298, 383 (1993); for detailed analyses and other references, see S. M. Bilenky, 
C. Giunti and W. Grimus, Prog. Part. Nucl. Phys. 43, 1 (1999); N. Okada and O. Yasuda, 
Int. J. Mod. Phys. A 12, 3669 (1997). 

[7] S. M. Bilenky, C. Giunti, W. Grimus and T. Schwetz, Phys. Rev. D60, 073007 (1999); B. 
Balatenkin, G. Fuller, J. Fetter and G. McGlaughlin, Phys. Rev. C 59, 2873 (1999); V. 
Barger, B. Kayser, J. Learned, T. Weiler, K. Whisnant, Phys. Lett. B489, 345 (2000); O. 
Peres and A. Y. Smirnov, [hep-ph / 001 1054| ; K. Balaji, A. Perez-Lorenzana and A. Smirnov, 



hep-ph/ 010 1005] R. N. Mohapatra, |hep-ph/0107274| ; Phys. Rev. D 64, 091301; K. S. Babu 



and R. N. Mohapatra, |hep-ph/020 17761 



[8] M. C. Gonzales-Garcia, M. Maltoni and C. Pena-Garay, |hep-ph/0108TJ7^ ; M. Maltoni, T. 
Schwetz, J.W.F. Valle, e-Print Archive: |hep-ph/0112"T03| . 

[9] M. Gell-Mann, P. Ramond and R. Slansky, in Supergravity, eds. P. van Niewenhuizen and 
D.Z. Freedman (North Holland 1979); T. Yanagida, in Proceedings of Workshop on Unified 
Theory and Baryon number in the Universe, eds. O. Sawada and A. Sugamoto (KEK 1979); 
R.N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[10] R. Foot and R. Volkas, Phys. Rev. D 52, 6595 (1995); Z. Berezhiani and R. N. Mohapatra, 
Phys. Rev. D 52, 6607 (1995); E. Ma and P. Roy, Phys. Rev. D 52, 4342 (1995); K. Benakli 
and A. Smirnov, Phys. Rev. Lett. 79, 4314 (1997); P. Langacker, |hep-ph/980528l| ; D. Sue- 



matsu, Phys. Lett. B392, 413 (1997); E. J. Chun, A. Joshipura and A. Smirnov, Phys. Rev. 
D 54, 4654 (1996); E. Ma, Phys. Lett. B 30, 286 (1996); Z. Chacko and R. N. Mohapatra, 
|hep-ph/9905388l 

[11] R.N. Mohapatra, Phys. Rev. D64 091301 (2001) 

[12] K.S. Babu, R.N. Mohapatra, Phys. Lett. B522 287 (2001) 



10 



[13] The Super-Kamiokande Collaboration, S. Fukuda et ai, Phys. Rev. Lett. 85, 3999 (2000). 

[14] A. Bandyopadhyay, S. Choubey, S. Goswami , K. Kar, Phys. Lett. B519, 83 (2001); G.L. 
Fogli, E. Lisi, D. Montanino, A. Palazzo, Phys. Rev. D64, 093007 (2001); J.N. Bahcall, M.C. 
Gonzalez-Garcia, C. Pana-Garay, JHEP 0108, 014 (2001); P.I. Krastev and A.Yu. Smirnov, 
e-Print Archive: [hep-ph/ 01081771 ; M.V. Garzelli and C. Giunti, JHEP 0112, 017 (2001). 

[15] S. Goswami and A. Joshipura, e-Print Archive: hep-ph/011027; P. Roy, S. K. Vempati, Phys. 
Rev. D65, 073011 (2002). 

[16] H.V. Klapdor-Kleingrothaus et a/., Mod. Phys. Lett. A37, 2409 (2001). 

[17] C. Marquet et ai, Nucl. Phys. Proc. Suppl. 87, 298 (2000). 

[18] E. Fiorini, Phys. Rep. 307, 309 (1998). 

[19] H.V. Klapdor-Kleingrothaus, J. Hellmig and M. Hirsch, J. Phys. G24, 483 (1998); L. Baudis 
et a/., Phys. Rep. bf 307, 301 (1998). 

[20] The MAINZ Collaboration, J. Bonn et ai, Nucl. Phys. Proc. Suppl. 91, 273 (2001). 

[21] The TROITSK Collaboration, V.M. Lobashev et ai, Nucl. Phys. Proc. Suppl. 91, 280 (2001). 

[22] The KATRIN Collaboration, A. Osipowicz et a/., |iep-ex/0109033| . 

[23] B. Achkar et ai, Nucl. Phys. B424, 503 (1995). 

[24] J. C. Pati and A. Salam, Phys. Rev. D10, 275 (1974); R. N. Mohapatra and J. C. Pati, Phys. 
Rev. D 11, 566 (1975); Phys. Rev. D 11, 2558 (1975); G. Senjanovic and R. N. Mohapatra, 
Phys. Rev. D 12, 1502 (1975). 



11 



